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ABSTRACT 

We present the first optical (590-890 nm) imaging polarimetry observations of the pre-transitional protoplan¬ 
etary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous 
studies. We detect the previously unseen far side of the disk gap, confirm the highly off-centered scattered-light 
gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from mil¬ 
limeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This 
new source of scattered light may contribute to the near-infrared interferometric signal attributed to the proto¬ 
planet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present 
a new model for the system architecture of LkCa 15 that ties these new findings together. These observations 
were taken during science verification of SPHERE ZIMPOL and demonstrate this facility’s performance for 
faint guide stars under adverse observing conditions. 

Keywords: circumstellar matter — planets and satellites: formation — protoplanetary disks — stars: individual 
(LkCa 15) — stars: pre-main sequence — techniques: high angular resolution 


1. INTRODUCTION 


Transitional disks are protoplanetary disks with heavily de 
pleted gaps or cavities i n their inner regions (e.g. Strom et al. 


1989 Calvet et aL||2005[ ). They are thought to represent a 
transitional state in disk evolution, in which the gas-rich pri¬ 
mordial disk grad ually disperses due to planet formation or 
othe r effe cts (e. g.|Bryden et al.|1999[|Owen & Jackson[2012[ 


[Alexander et al.l20l4|). 

One transitional disk target that has garnered particular at¬ 
tention in recent years is LkCa 15, a 3-5 Myr old ~1 M© 
star located approximately 140 pc away ( Simon et al.} 2000| 
in the Taurus-Auriga region. A wide gap in the disk extend¬ 
ing out to ~50 AU can be inferred from the infrared s pectral 
energy distribution (SED; |Espaillat et al.|2007[|2008| , but it 
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has also been spatially r esolved both in scatter ed light at near- 
infrared wavelengths ( [Thalmann et al.l|2010|) an d in thermal 
emission at mm/sub-mm wavelengt hs fPietu et al.|2007| An 


.|2Q 11 1 Jlsella et al.|2014 ). In near-infrared images, 
Ige of the gap appears off-center from the star, sug- 


drews et al. 

the outer edge of 1 _ __ 

gesting an eccentric gap ([Thalmann et al.|20>T4] ), although the 


mm/sub-mm images are consisten t with a circ ular gap ([An 
drews et al.|2011[|Isella et al.|2014| ). In addition, [Kraus & Ire¬ 


land ^2012]rdiscovered a protoplanet candidate at ~ 16-21 AU 


using Sparse Aperture Masking (SAM). While some follow 
up observations onl y provided upper limits ( Isella et al.|2014 


Whelan et al.|2015 ), newest results from H a differential imag 
l th( 


ing confirm the accreting protoplanet hypothesis (K. Follette, 


p.c.). 

Here we present new scattered-light images of the LkCa 15 
disk taken at visible wavelengths using Polarimetric Differ¬ 
ential Imaging (PDI). The observations provide better spatial 
resolution and a more robust determination of the distribution 
of reflected light from the disk than previous data based on 
infrared Angular Differential Imaging (ADI), in which for¬ 
ward modeling is necessary to account for self-subtraction ef- 
fects ([Thalmann et al.|201 l[|Milli et al.|2012[|Thalmann et al.| 
[2013]). 


2. OBSERVATIONS AND DATA REDUCTION 


LkCa 15 was observed on two nights during Science Verifi¬ 
cation Time for the newly commissioned high-contrast imager 
SPHERE (Spectro-Polarimetric High-contrast Exoplanet RE- 
search; Beuzit et al. |2008| ) at the European Southern Obser¬ 
vatory’s Very Large Telescope facility. Its visibl e-light instru- 


al. 2008 [ Schmid etal. 2012 

), was used with the extreme adap- 

tive optics system (SAXO; 

Fusco et al. 2006) to obtain deep 


the R and I bands (590-890 nm). ZIMPOL offers extreme 
polarimetric contrast ( [Schmid et al.|2012| at a pixel scale of 
3.6 mas x 7.2 mas and a field of view of 375 x 375. 

The first observation, on 2015-02-02, comprised seven po¬ 
larimetric cycles of four half-wave plate positions (Stokes +Q , 
-Q, +£/, -£/), each exposed for 8 x 20 seconds at a time, for 
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a total integration time of 75 minutes. The second observation 
was executed in the same way on 2015-02-12. 

The exposure time had been chosen to saturate the star’s 
point-spread function (PSF) to maximize sensitivity at larger 
separations. However, the performance of the adaptive op¬ 
tics was poor and highly variable throughout both observa¬ 
tions, which rendered the star unsaturated in most frames and 
yielded a full width at half maximum (FWHM) of ~35 mas 
(« 2A/D). This can be ascribed to LkCa 15’s apparent mag¬ 
nitude of R « 12 mag at the faint limit of SAXO’s guide 
star specifications, combined with unfavorable weather con¬ 
ditions and high airmass (~1.5). Furthermore, a “comet-tail” 
aberration of the PSF core is observed in some frames, likely 
caused by malfunction of the atmospheric dispersion correc¬ 
tor (ADC). 

After correcting the raw ZIMPOL data for cosmetics and 
field distortioicwe applied the PDI data reduction method de¬ 
tailed in Avenhau s et al.| ( |2014a| ). It includes the following 
steps for reducing instrumental polarization: (1) flux equal¬ 
ization of each simultaneous pair of raw images in order to 
remove the polarization of the stellar PSF, (2) double-quotient 
calculation of Stokes Q and U polarized flux images for each 
polarimetric cycle, (3) empirical correction of polarimetric 
throughput differences in Stokes Q and U , and (4) empirical 
correction of angular misalignment of the half-wave plate. 

Lastly, we transform the Cartesian coordinate system of 
the Stokes formalism ( Q , U) into a polar coordinate system 
(Q^Uf) in which positive Q$ polarization is defined as az¬ 
imuthal with respect to the star ( [Benisty et al.||2015| ). Under 
the single-scattering assumption, the scattering polarization 
should appear in the Q$ image (typically as positive signal), 
whereas the U$ image should contain no scattered-light sig¬ 
nal a nd can be used to estim ate the noise level in the observa¬ 
tion ( [Avenhaus et al.|2014a| ). While this assumption does not 
strictly hold for inclined, optically thick disks, as is the case 
for LkCa 15, the resulting deviations are small and automat¬ 


ically accounted for in the error estimation (Avenhaus et al. 
12014b]) . 


Since we find no significant differences between the two 
epochs, we combine both data sets into one final set of Stokes 
(/, Qf, U<p) images, which improves the quality of the results. 
Since the data are sensitivity-limited, we apply no frame se¬ 
lection. We convolve the polarization images with a circu¬ 
lar aperture of 7 pixels (= 25 mas) in diameter to reduce shot 
noise. 


3. RESULTS 
3.1. Imagery 

Our RI -band imaging polarimetry of LkCa 15 is shown in 
Fig. [I] Despite the bad PSF quality, the instantaneous nature 
of dual-beam polarimetry allowed us to remove the starlight 
effectively and resolve the LkCa 15 disk at high contrast and 
unprecedented angular resolution. 

The Qtf, image reveals a crescent of polarized flux trac¬ 
ing the shape of the near-side edge of the o uter dis k seen 
in infrared angular ADI observations ( [Thalmann et al.|[2010[ 
|2014| , and a previously unknown bright area of polarized flux 
surrounding the star at small separations (< 0725, see panel 
[lj)). Both features appear consistently positive in the Q$ im¬ 
age and are significantly brighter than the residuals in the 
image at similar radii, which identifies them as scattered light 
from circumstellar material. 

Fig. [I]: provides signal-to-noise (S/N) maps produced by 





Figure 1 . Imaging polarimetry of LkCa 15 based on SPHERE ZIMPOL ob¬ 
servations at 590-890 nm. In each panel, the left and right sides show the Qq 
and polarized flux images, respectively. Assuming single scattering, the 
disk appears as positive flux in the Q$ image, while the U$ image provides a 
noise estimate. All images are convolved with a 7-pixel (=25 mas) diameter 
circular aperture to reduce shot noise. The small insets show the shape of the 
PSF core. The location of the star is marked with black-and-white circles. 
The position angle of the North direction is 30°. The black triangles indicate 
the zero level in the color scales, (a) Field of view of 178 at a linear stretch 
with a maximum of ~10 -4 times the brightest pixel of the stellar PSF (which 
is lightly saturated in some of the exposures). The known crescent of scat¬ 
tered light from LkCa 15’s outer disk is visible, along with a newly discovered 
inner disk component close to the star, saturating the color scale. The white 
corners delimit the area used for (b). (b) The central 07 6 shown at 3x magni¬ 
fication at a linear stretch with a maximum of ~10 -3 times the brightest pixel 
of the stellar PSF, showing an unsaturated view of the inner disk component. 
The cros shairs indicate the reported K -band location of LkCa 15 b in [Kraus] 
| & Ireland) {201 2| . The dashed circle marks the noise-dominated central area, 
(c) Signal-to-noise (S/N) maps at a stretch of [-5,5]cr. The same noise map 
is used for both the Q$ and Uq images. Both the outer disk and the new inner 
disk component are detected at high significance, (d) The 178 field of view 
with each pixel multiplied by the square of its de-projected distance from the 
star, assuming an inclination of 50°. 
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Table 1 

Numerical results. 


Figure 2. Ellipse fits to the LkCa 15 disk gap edge. The underlying image 
is the r 2 -scaled Q$ image as in Fig. [ill displayed at a harder stretch, whereas 
the light blue ellipse represents ournt to its gap (see text). The location of 
the star is marked with a star symbol, (a) Visual confirmation of fit for the 
blue ellipse, (b) Modified maximum merit fit to the best-fit model from the 
ADI-based forward-modeling analysis injThalmann et al.||2014) overplotted 
in yellow. The ellipse is registered against the underlying image by its star 
position, (c) Fit ellipse (green) to the gap edge in 7 mm emission as reported 
in Fig. 1 of|Isella et al.|(|2014|, registered by the assumed star position, (d) 
The same ellipse as in (c) aligned with our blue fit ellipse, demonstrating that 
the shapes match well. The small green circle marks the true position of the 
star relative to the green ellipse. 

dividing the pixel values of star-centered annuli in the Q$ and 
Ucj) images by the standard deviation of the U$ image in the 
same annuli. Both features achieve S/N values above 5<x over 
a large number of resolution elements and are therefore de¬ 
tected at high confidence. The inner disk is detected down to 
a separation of ~1 FWHM (36 mas). 

To estimate the dust distribution from the scattered-light in¬ 
tensity, we compensate for the r ~ 2 fall-off of the illumination 
of dust particles with physical distance r from the star. We cal¬ 
culate r for each pixel by deprojecting its apparent distance 


Fit ellipse to the disk gap edge 

in 




in mas 

in proj. 

AU 

Semimajor axis a 

338 

(+11,-18) 

47.4 

(+1.5, -2.5) 

Semiminor axis b 

245 

(+21,-11) 

34.3 

(+3.0,-1.5) 

Major-axis offset x 

-3 

(+11, -14) 

-0.5 

(+1.5, -2.0) 

Minor-axis offset y 

-79 

(+11,-18) 

-11.1 

(+1.5, -2.5) 

Characteristic length scales of the inner scattering component 


in mas 

in proj. 

AU 

Major axis (East) /e 

69 

(+15,-11) 

9.7 

(+2.1,-1.5) 

Major axis (West) /w 

74 

(+11, -9) 

10.4 

(+1.6,-1.2) 

Minor axis (North) /n 

51 

( +6, -5) 

7.2 

(+0.9, -0.7) 

Minor axis (South) 1$ 

38 

(+10, -6) 

5.3 

(+1.6, -0.9) 


et al. 12012} |Thalmann et al.||2013| ) and ther efore require for¬ 
ward modeling to interpret their results (jThalmann et al. 


mg to interpret their results ( jThalmann 
2014| ). Imaging polarimetry, on the other hand, yields mor¬ 


phologically sound rep resentations of the disk ’ s appearance 
in polarized light (e.g., Hashim oto et al.||2011} Quanz et al. 


20i3a|b[|Garufi et al.|2013[|Benisty et al.|2015|i. 

We characterize the sharp gap edge in the r 2 - scaled Q$ im¬ 
age by fitting it with an ellipse using a modified versio n of 
the m aximum merit method introduced in IThalmann et al.l 
( 2011| ). We generate a number of ellipses with semimajor 
and -minor axes a , b displaced from the star’s position by off¬ 
sets x, y along the two axes, respectively. For a parameter 
set (< a , b , x,y), we consider the two concentric elliptical annuli 
bounded by the triplet of ellipses (a-S,b-^S,x,y),(a,b,x,y), 
and (a + S, b + x, y) with an annulus width of S = 10 pix¬ 
els (= 36 mas = 1 FWHM). The merit function rj(a , b , x,y) is 
defined as the mean pixel value in the outer annulus of the r 2 - 
scaled Q$ image minus the mean value of the inner annulus. 
The best fit is achieved when the merit function is maximized, 
i.e., where the mean brightness increase from the inner to the 
outer annulus is greatest. As a measure of uncertainty, we 
define the “well-fitting” family of ellipses by a threshold of 
tj > 7/ max - CTjj, where cr v is the standard deviation of merits of 
ellipses in the U# image. 

As demonstrated in Fig. [2^, the best-fit ellipse matches the 
visual locus of the gap edge convincingly. Its numerical pa¬ 
rameters and uncertainties are provided in Table [I] Figure (2J) 
compares this fit to a second ellipse obtained by applying 
the maximum merit method to the best-fit model from our 



dence of scattering polarization on the phase angle. However, 
we note that phase angles should change by at most ~6° along 
radial paths on the visible disk surface; thus, the location of 
the gap edge should remain largely unaffected. 

The eastern ansa appears brighter than the western one on 
the 1 cr level (evaluated in the r - scaled image in 50-pixel di¬ 
ameter apertures). 

3.2. Geometric analysis 

ADI observations of circumstellar disks suffer from flux 
loss and oversubtraction effects ( jThalmann et al.|20TT} |Milli| 


significant mismatch with our off-centered ellipse. While the 
two ellipses can be brought to coincide well with each other 
by translating one of them by y « 80 mas (Fig.|2]l), this offset 
is highly significant compared to the astrometric accuracy of 
the 7 mm interferometry (10mas,|Isella et al.|2014j) and of the 
ZIMPOL data (< 10 mas) and must therefore be physical. 

3.3. Surface brightness profiles 

Figures [3t and[3]3 plot the surface brightness profiles in the 
Q(P image (without r 2 scaling) along the major and the mi¬ 
nor axis, respectively. The plotted values are pixel values as a 















































4 


Thalmann et al. 



-0.5 0.0 0.5 

Separation along major axis (arcsec) 



-0.5 0.0 0.5 

Separation along minor axis (arcsec) 


Figure 3. Polarized surface brightness profiles for the LkCa 15 disk, (a) 
Profile of the Q ( p image along the major axis averaged over a width of 21 
pixels (=76 mas) centered on the star (thick black curve). The family of thin 
grey curves are obtained by adding similar profiles taken in the U# image at 
various position angles to the black curve, and therefore visualize the noise 
level in the data. The contrast of a given pixel is calculated by the ratio of its 
value and the maximum pixel in the stellar full-intensity PSF. The normaliza¬ 
tion is approximate due to the stellar PSF being saturated in some exposures. 
The blue dashed lines mark the location of the light blue fit ellipse in Fig. [2] 
The red lines are exponential fits to the inner scattering component. Positive 
separation values correspond to the western side of the disk. The size of the 
PSF’s FWHM is given for comparison, (b) Profile of the Q image along the 
minor axis averaged over a width of 81 pixels (= 292 mas) centered on the 
star. Positive separation values correspond to the northern side of the disk. 


function of separation averaged over a width of 21 pixels (= 
76 mas) and 81 pixels (= 292 mas), respectively, and normal¬ 
ized by the brightest pixel value of the star’s intensity PSF. 
Since some exposures in the data are saturated, the normal¬ 
ization is only approximate. 

The new inner disk component appears as a roughly expo¬ 
nential rise of surface brightness within a radius of ~0725. 
We fit each side of the structure in each of the two profiles as 
proportional to exp (-6/1), with 6 the angular separation from 
the star and / a characteristic length scale, using least-squares 
fitting. These fits are listed in Tableland illustrated in Fig. [3] 

To gauge the uncertainty of these results, we extract a num¬ 
ber of surface brightness profiles from the noise-like U# im¬ 
age at various position angles, add each of them to the science 
profiles from the Q$ image, and recalculate the exponential 
fits. The standard deviations of these “disturbed” l values are 
then adopted as the error bars for the table values of /. The dis¬ 
turbed surface brightness profiles are shown in grey in Fig. [3] 

The polarized intensity is likely underestimated at the in¬ 


nermost separat ions due to overlapping PSF wings in the Q 
and U images ( [Avenhaus et al.|2014b| ). 


4. DISCUSSION 
4.1. Inner disk component 

The newly discovered inner disk component, detected from 
~30 AU down to the inner limit of ~7 AU (Fig. [3]), overlaps 
with the location of the protopla net candidate at -16-2 1 AU 
derived from SAM observations ( |Kraus & Ireland|2012| ). As 
with the outer disk, we expect anisotropic scattering to pro¬ 
duce a distinct full-intensity forward-scattering maximum on 
the near side, which coincides with the position angle of 
the protoplanet candidate. Taken together with the forward¬ 


scattering maximum of the outer disk edge (Thalmann et al. 
2014| ), these disk features likely interfere with the reconstruc 


tion of SAM data, as in the case of T Cha ( [Olofsson et al.| 


2013| ). However, given the recent H a detection (K. Follette, 


p.c.), a proto planet embedded in the inner disk similar to 
HD100546 b ( [Quanz et al.||2Q13a) ) remains the most plausi¬ 
ble interpretation. 

The inner disk component also sheds light on the extent of 
the gap as inferred from previous non-resolved observations. 
[Espaillat et al. ( 2007| proposed a three-component structure 
with an outer disk at ~50 AU consistent with the millimeter 
cavity, an inner disk starting at the dust sublimation radius to 
explain the near-infrared flux, and an optically thin compo¬ 
nent out to ~5 AU to explain the prominent silicate emission 
feature. Our detection of scattered light between 7-30 AU 
shows illuminated dust is present at these locations above the 
disk mid plane. An extrapolation of the gap wall temperature 
(T « 100 K/ ^a/a wa n) yields a temperature of 150-250 K, 
consistent with the presence of the 10 and 20 pm silicate emis¬ 
sion feature in the SED. We hypothesize that this new scat¬ 
tering structure represents the optically thin part of the inner 
disk extending out to at least 30 AU. A quick adaptation of 
the disk model from [Mulders et al.[ ( |2010[ ) indicates that this 
can be reconciled with the SED, since the contribution of the 
optically thin part to the SED is minor. This geometry (nar¬ 
row gap, depleted inner disk, millimeter cavity) qualitatively 
matches the imprint of a single planet on a disk (e.g. |Zhu et| 
al.|2012HJang-Condell & Turner|2012) . 


4.2. Off-centered gap 

Our RI -band polarimetry confirms the prediction from K- 
band ADI forward-modeling ( Thalmann et al.|2014| ) that the 
outer boundary of the gap in the LkCa 15 disk appears off- 
centered in reflected light, with a highly significant offset of 
y ~ 80 mas (a third of the gap’s projected radius) between 
the star and the gap’s apparent center. While disk gaps are 
known to vary in size as a function of observed wav ele ngth in 


systems with planet-induced dust filtering (de Juan Ovelaret 
al.|2013 ), a lateral displacement of the gap is more difficult to 


justify and has not been predicted by numerical simulations. 

The outer disk’s vertical structure cannot account for this 
displacement since the gap wall profile is tapered (‘round’) 
rather than vertical ( [Thalmann et al.|[2014| ). This is qualita¬ 
tively confirmed by the wide, homogeneous distribution of 
scattered light in Fig.[l|i 

Another explanation is provided by shadowing. |Espaillat et| 
originally proposed an inner disk coplanaFwith the 
c, casting a s hadow onto the outer disk’s inner wall 


aL|([2007 

outer 3is 


Muld ers 
et al. 


ders et al7| ( |2010j ). Our forward modeling in Thalmann 
.|( |2014 ) achieved significantly better fits with an unshad- 
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Figure 4. Schematic of the proposed architecture of the LkCa 15 disk (not 
to scale). Red lining marks the disk surfaces illuminated by the star. Dotted 
black lines trace the edges of the shadow cast by the inner disk. 


owed gap wall, which led us to speculate that the inner disk— 
and its shadow—were tilted out of the outer disk’s plane. On 
the other hand, this arrangement should result in dark lanes 
cutting across the outer disk where the shadow plane inter¬ 
sects it, as seen in HD 142527 ( |Marino et al.|2015| ). No such 
lanes are apparent in the LkCa 15 imagery. 

This conundrum can be resolved if the inner disk has a 
slightly higher inclination than the outer disk (Fig. 0]). The 
near-side disk rim as observed from Earth would then "be fully 
exposed to stellar irradiation, whereas the far-side disk sur¬ 
face would be obscured at the rim and only become visible 
further out where the disk surface flares out of the shadow. 
This would create an off-centered apparent gap in reflected 
light even wit h a symmetric physical disk gap. For the disk 
architecture in Thalma nn et al.| ( 20l4] ), the requisite differen¬ 
tial tilt would be of order 1°, which could be i nduced by a 
planet on a tilted orbit (cf. Mouillet et al.|1997[|Augereau et| 


al.|2001[[Lagrange et al.|201^~ 

A similar effect could be caused by a shadowing structure 
vertically extending from an otherwise coplanar inner disk 
( [Espaillat et al.|2011| , though this would imply dramatic vari¬ 
ability on the timescale from days to years. While we cannot 
exclude this possibility^we find no significant variability in 
our observations ( [Thalmann et al.[2014t this work). 
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